lifecontingencies (version 1.3.1)

multiple life probabilities: Functions to deals with multiple life models

Description

These functions evaluate multiple life survival probabilities, either for joint or last life status. Arbitrary life probabilities can be generated as well as random samples of lifes.

Usage

exyzt(tablesList, x, t = Inf, status = "joint",  type = "Kx", ...)
pxyzt(tablesList, x, t, status = "joint", fractional=rep("linear",
length(tablesList)), ...)
qxyzt(tablesList, x, t, status = "joint", 
fractional=rep("linear",length(tablesList)), ...)
rLifexyz(n, tablesList, x, k = 1, type = "Tx")

Arguments

tablesList

A list whose elements are either lifetable or actuarialtable class objects.

x

A vector of the same size of tableList that contains the initial ages.

t

The duration.

n

The size of sampled life duration matrix.

status

Either "joint" (for the joint-life status model) or "last".

type

"Tx" for continuous, "Kx" for curtate.

fractional

Fractional lives assumption.

…

Options to be passed to pxt.

k

Fractional frequency option.

Value

An estimate of survival / death probability or expected lifetime, or a matrix of ages.

Details

These functions extends pxyt family to an arbitrary number of life contingencies.

References

Broverman, S.A., Mathematics of Investment and Credit (Fourth Edition), 2008, ACTEX Publications.

See Also

pxt,exn

Examples

# NOT RUN {
#assessment of curtate expectation of future lifetime of the joint-life status
#generate a sample of lifes
data(soaLt)
soa08Act=with(soaLt, new("actuarialtable",interest=0.06,x=x,lx=Ix,name="SOA2008"))
tables=list(males=soa08Act, females=soa08Act)
xVec=c(60,65)
test=rLifexyz(n=50000, tablesList = tables,x=xVec,type="Kx")
#check first survival status
t.test(x=apply(test,1,"min"),mu=exyzt(tablesList=tables, x=xVec,status="joint"))
#check last survival status
t.test(x=apply(test,1,"max"),mu=exyzt(tablesList=tables, x=xVec,status="last"))
# }